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Abstract. This is a mainly expository sketch showing how some integrable systems (e.j 
KP or KdV) can be viewed as quantum mechanical in nature. 



1. INTRODUCTION 

We review briefly some ideas and constructions for dKP (from KP) and its Hamiltonian 
dynamics, discuss its quantization via Moyal back to KP, and treat the KP hierarchy itself as 
a sequence of quantum Hamiltonian systems based on the dKP dynamics. Some connections 
to q-theories are also indicated. The constructions are all well known (cf. [HI IH1 ini El El 

rra rrrn run rm rrni nui f^i f^m n^i r?rn r?Hi fm r^i r^i rrrn isa esd but we try to organize 

matters and perhaps see them from a different perspective. We will use dKdV and KdV as 
a source of simple examples. 

2. BACKGROUND 

There is an enormous literature on KP and KdV which we do not try to reference (see e.g. 
[Ill 1221 |2*§] ). For a brief sketch we follow |17j and simply write down the relevant formulas. 
Thus one begins with two pseudodifferential operators (Al) L = d + ^i° n n+i9 _n and 
W = 1 + Y1T w n d~ n (Lax and gauge operator respectively) with L = \VdW~ 1 . Note here 
the generalized Leibniz rule 

(2.1) d n f = jr( l \d k f)d n - k 

^ ' 

The KP hierarchy is then determined by Lax equations (A2) d n L = [B n ,L] = B n L — 
LB n where d n ~ d/dt n and B n = 17] is the differential part of L n = U\ + L™. One 
can also express this via the Sato equation (A3) d n WW~ l = — L™. Now define wave 
functions via (A4) ip = Wexp{£) = w{t,X)exp{£) where £ = X^i n A n an d w(t,X) = 
1 + w n (t)\~ n (with t\ = x). There is also an adjoint wavefunction (A5) ip* = 
W*~ l exp(-£) = w*(t,\)exp(—£) where w*(t,X) = 1 + (£)A -< ; there are then equa- 

tions 

(2.2) Lip = XtP; d n ip = B n iP; L*iP = Xip*; d n ^* = -B*^* 
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Next one has the fundamental object, namely the tau function, which yields the wave 
functions via vertex operators X, X* in the form 

(2 - 3) A) = —^T = — W) — ' * ( 5 } = W) = W) 

We write r±(t) = r(t ± [A" 1 ]) = exp(±£(<9, A _1 )t where d = (d u (l/2)<9 2 , (l/3)<9 3 , • • • ) and 
t ± [A -1 ] = (ti ± A _1 ,i2 ± (1/2)A -2 , • • • ) and we can also specify 

(oo \ oo 

Y.tnX 11 ) =Y J P J (h,--- ,tj)\ j ; P j(t) -- 

where the pj are elementary Schur polynomials (note Ylj^j = n )- O ne recalls also the 
famous Hirota bilinear identity (A6) § ip(t, X)ip*(t' , A)d\ = (residue integral around oo). 
Using (JT3J) this can be written as (A7) ^ r(t-[A~ 1 ])r(t'+[A~ 1 ])ea:p[^(t, A)-^, A)]dA = 
which leads to the characterization of the tau function via (t — > t — y and i ; — ► t' + y) 

(2.5) ^ Pn (-2yK +1 (D)exp ^f>AU t • r = 

where A is the Hirota derivative defined as (A8) D^a-b = (d m /ds™)a(tj + Sj)b(tj — Sj)| s =o 
and D = (D\, (1/2)D 2 , (1/3) A, • • • )• ^ n particular from the coefficient of the free parameter 
y n in (|2.5|) one obtains (A9) D\D n T ■ r = 2p n+ i(D)r ■ r and these are called Hirota bilinear 
equations. Motivated via finite zone situations and Riemann surfaces (cf. |22j 'l where 
the tau function is intimately related to theta functions one can express the Fay trisecant 
identity in a form referred to as the Fay identity (A10) ^2 cp (so — s%)(s 2 — S3)r(t + [sq] + 
[si])r(i + [^2] + [S3]) = where c.p. means cyclic permutations (cf. ^[^2j). Differentiating 
this in so, setting so = S3 = 0, dividing by s\s 2 , and shifting t — > t — [S2], leads to the 
differential Fay identity (d = d x ) 

(2.6) T(t)dr(t + H - [s 2 ]) - r(t + [ Sl ] - [s 2 ])dr{t) = 

= (sr 1 - a?){r(t + N - [s 2 ])r(t) - r(t + [ Sl ])r(t - [s 2 ])} 
The Hirota equations (A9) can also be derived from l|2.6|) by taking the limit s\ — > s 2 . 

For the dispersionless theory (dKP) one thinks of fast and slow variables for example 
taking t n — > et n = T n and t\ = x — > ex = X. Then letting e — > the KP equation 
(All) u t = (l/4)u xxa; + 3uu x + (3/4)0- 1 nj / j / (u ~ u 2 ) goes to d T U = 2UU x + {3/4)d- 1 U Y Y 
where d ~ dx now (for a discussion of the passage u n — > J7 n we refer to ^3 - generally 
we think here of u n (T/e) ~ U n (T) + 0(e) etc.). Now take a WKB form for the wavefunction 
of the form (A12) ip = exp(S(T, A)/e) and define P = dxS 1 ; then as e — > the equation 
Lip = AV> becomes (A13) A = P+£~ U n +iP~ n with inverse P = Pi+l^- Further 

from B n *P = YA) b nm{ed) m ^ one obtains (A14) d n S = <B n (P) = A^ where B n = L% -> 
25n = A" = X^o bnmP m - Consequently the KP hierarchy becomes (A15) <9 n P = <923 n (note 
9 n S = ^B n =>■ 9 n P = 953 n and 3 n ~ d/dT n here). One can now write 

°° °° r) TP / 1 \ 

(2.7) 5 = £ T„A« - ^ ~~^~~ m > r = exp (L F (T) J 
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and there results P = dS = 58 1 with 



i 

oo 



(2.8) <B„ = d n S = X n - V ^H\~ m 

£■ — ' m. 



m 
l 



Next following I5n| one can manipulate the differential Fay identity in various interest- 
ing ways and in [TJ\ the author and Y. Kodama derived a dispersionless limit of the Hi- 
rota equations which provided polynomial identities among the coefficients F mn = d n d m F. 
These dispersionless Hirota equations can be written in the form (A16) Fij = pj + i(Z\ = 
Q,Zz,--- ,Zj + %) where p n ~ Schur polynomial and Zj = Yl m +n=j(^mn/mn). One shows 
also that the F mn can be expressed as polynomials in the Pj+i = F\j/j. There are many 
other results in (including connections to D-bar techniques) and dKP has been de- 
veloped and used in a number of fascinating contexts in recent years (we refer e.g. to 

REMARK 2.1. Another important aspect of dKP is the Hamiltonian theory origi- 
nating in (cf . also [111 1151 1161 117j ) . It is convenient here to rescale the variables via 
t n — > T' n = nT n with Q n = {l/n)% n . There results 

(2-9) P' n = = dQ n ; X' n = — = -d P £} n 

One can also show that (cf. |17j ) 

(2.10) PM _ P(A ) °£^1%- 

which is in fact equivalent to the dispersionless differential Fay identity 

(2.11) ^-mx-n^rnn = { L _ y fM~ n - X~ n F^A 
m,n=l x ' 

The kernel in (|2.10|) represents a Cauchy type kernel and following Kodama has a version 
on Riemann surfaces related to the prime form. ■ 

3. INTEGRABLE SYSTEMS AND MOYAL 

For background we follow |11[ 1121 113j and one recalls for wave functions tp there are 
Wigner functions (WF) given via 

(3.D /(.,,) -i/W.-^W-W. + 5» 



Then defining f * g via 

(3.2) f*g = f eX p 



9- 



f(.r.p) * <)(-r. pi ~ f j X + ^d p ,p- y d x ) <](.r.p) 
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time dependence of WF's is given by (H ~ Hamiltonian) 

(3.3) d t f(x,p,t) = -rz{H * f(x,p, t) - f(x,p,t) * H) = {H,f} M 

in 

where {/, g}M ~ Moyal bracket. As H — > this reduces to dtf — {H, /} = (standard 
Poisson bracket). One can generalize and write out (J3.2|) in various ways. For example 
replacing ih/2 by k one obtains as in |36| 

(3-4) / * g = E ^ ( ■ )(9id s p - j fm- j di9) 

leading to ({/, g} K = (/ * g - g * /)/2«) 

,2s 2.+1 /2s + 1 



(3-5) {/, g}, = £ ^1)! £ ( / J C^ 1 ^/)^ 



o ^ ' j=0 

which can also be utilized in the form 



(3.6) / * g = fe r(0*%-*ap3.) g = e^ d ^- d ^ d ^f(x 1 ,p 1 )g(x 2 ,p 2 )\^ p) = 
_°° ( lY^r+s ar+s f fvr+s n K n( -\\n— s 

= y i 1J " ° J ° J. = y lj (d n x - s d s j) (d s x d;- s g ) 



(| r\s\ dx r dp s dp r dx s ^—^ s\(n — s)\ 



oo „ n 
K 



E^- £(-*)' W" r /) (srqM 



n 





Note e.g. 

(3.7) g* f = g(x + Kd p ,p- nd x )f = f(x - nd p ,p + nd x )g 
The Moyal bracket can then be defined via 

(3.8) {/, g} M = -{/ Sin[nCd x ~d p - *d p ~d x )]g} = — (/ * g - g * /) = 

= E^t9 E(-if ( "7 1 ) >-*r' J .n>;r^ ^: 

corresponding to k — ► in in (|3.12|) below. 

Thus in 36 (cf. also [53]) one writes the Sato KP hierarchy via (v-2 = 1, w_ i = 0) 

oo 

(3.9) d m L = [L™, L] (m > 1); X = £ ^(x)^-"" 1 

for x = (x, ■ ■ ■ ) while the Moyal KP hierarchy is written via (u-2 = 1, u_i = 0) 

oo 

(3.10) A ~ A = y j u n {x)p- n - 1 - d m A = {A™, A} M (m > 1) 
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where A™ ~ (A* m ) + with 



(3-11) / * 9 = E ^ E(-!) j 

s - i=o VJ/ 

leading to 

00 2s 2s+l /„ , 1 \ 

(3.12) {/, g} K = E T^VW E ( mtip+^fl^-Wpg) 

+ ^- i=0 \ J / 

Note lim K ^ {f, g} K = {f,g} = f x g x - f x g\ so {KP) M difP as « -» 0, namely <9 m A = 
{A™, A} with A m ~ A---A. The isomorphism between (KP)sato and (KP)m is then 
determined by relating i> n and it n in the form (k = 1/2) 



(3.13) n n = E 2 " 



3 n ~ 3 



where n = 0, 1, • • ■ and u J = c^vo ( see Remark 3.1 for enhancement). 

REMARK 3.1. This can be further clarified as follows. First we recall an important 
paper [HHj where one considers star products of the form (Bl) f*g = fg + Yln>i h n B n (f, g) 
with bilinear differential operators B n . In particular in [3S] one shows that any bracket of 
the form 

00 00 r s 

(3.14) {f,g} = EE Ar+ " 2 EE 6 ^(^^/)(^^) 

r=l s=l j=0 k=0 

may be transformed to one with 600,10 = 600,11 = and any such bracket satisfying the 
Jacobi identity must be of the form 

00 r s 

(3.15) {f,g} ^r^^lri*^- 3 '/)^) 

r=l j=0 k=0 

By suitable calculation one shows also that (|3.15|) plus Jacobi is equivalent to Moyal. Note 
that the Jacobi condition for {/, g} = (l/h)(f * g — g * /) can be proved directly via 
associativity of * (exercise). Thus (B2) {{/, g},h} + {{h,f},g} + {{g,h},f} = 0. Now 
to connect the dKP theory with bracket (|3.12|) to a k — KP theory with PSDO bracket 
consider the PSDO symbol bracket (B3) AoB = £(l/fc!)<9|.4(x, £)d*B{x, £) (cf. QUEE1) 
where A ~ ^ a^x)^, d^A = ]T aj(x)<9|£\ and d*A = ^x a i{ x )C- Note also 

(3.16) Ao K B = Ae Kd ^ 9x B = Y,^jd^Ad^B 
and the bracket based on this; thus 

(3.17) A o K B = A(x, £ + Kd x )B{x, £); B o K A = B(x, £ + K d x )A(x, £) 

and (B4) {1/k)A o K B - B o K A) = {A, B} K is of the form (l3~TBl with b rr0 / 0, b r0r + 0, 
and all other coefficients equal 0. Also &no = —6101 and the Jacobi identity will follow 
from associativity so in fact any bracket such as (B4) is equivalent to Moyal in the symbols 
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involved (k is arbitrary). Note here that associativity is not obvious but is proved in |48j 
(note the k can be absorbed in £ by rescaling). The trick is to use the formula 

(3.18) (aO(ba = «E ( n ^)d k M n ~ k+r 

which shows that (for k = 1) 

(3-i9) ^ = E^c^ = £^- 

• £ *nnin - 1) • ■ ■ {n - m + ■ £ bf 1 ^ = £ a n Q 6<™ 

Now for associativity one checks that (B5) [£ n (a£ r )]fr = £ n [(a£ n )&] and we refer to [IS] f° r 
further details. ■ 

We recall also the standard symbol calculus for PSDO following e.g. |36l 1501 I58j (cf. 
(B3), (j3,16|) . (|3.17|) ). First one recalls from [SEj the ring or algebra 21 of pseudodifferential 
operators (PSDO) via PSD symbols (cf. also [SH] for a more mathematical discussion). 
Thus one looks at formal series (B6) A(x,£) = J^-oo a i( x )C where £ is the symbol for d x 
and di(x) € C°° (say on the line or circle). The multiplication law is given via the Leibnitz 

rule for symbols (B7) A(x,£) o B(x,£) = Efc^AO^f B ^ where = 

E n -oo a i( x )(^) {k) and B { x k \x,£) = E-oo b i k H x )? with b i k) W = d k h(x). This gives a Lie 
algebra structure on 21 via (B8) [A, B] = A o B — B o A. Now let A be a first order formal 
PSDO of the form (B9) A = d x + Oi(jc)Sj where £ ~ (s, i 2 , t 3 , ■ • • ). Then the KP 

hierarchy can be written in the form (BIO) (dA/dt m ) = [(A m ) + , A] which is equivalent to 
a system of evolution equations (Bll) (da,i/dt m ) = fi where the fi are certain universal 
differential polyomials in the ctj, homogeneous of weight m + \i\ + 1 where a?^ has weight 

|i| + l for a? ~ d x a. Somewhat more traditionally (following - modulo notation and 
various necessary analytical details), one can write 

(3.20) Au{x) = {2-k)- 1 j e ix< a(x,£)u(£)dZ 

where = J exp(—ix ■ £)u(x)dx. One takes D = (l/i)d x and writes a = symb(A) with 
A = op(a) ~ A where the • is to mod out (we will not be fussy about this and will 

simply use A) . The symbol of A o B is then formally 

(3.21) (a®b)(x,0 =J2^<^0D^b(x,0 

corresponding to (B7), while [A, B] = AB — BA corresponds to the symbol (B12) {a, b} = 

(da/d()(db/dx) - (da/dx) (db/d£) (note P{D)T = P(£)f). In any event it is clear that 
the algebra of differential operators on a manifold M (which we have sometimes loosely 
referred to as quantum operators) may be considered as a noncommutative deformation of 
the algebra of functions on T*M and the extension to PSDO brings one into the arena of 
integrable systems. 
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4. QUANTUM MECHANICS 

We refer now to jHS] where a lovely survey appears (cf. also ^2 EE1 HH EEl E2] ) • There 
are three logically autonomous alternative paths to quantization, namely Hilbert space 
operators, path integrals, and deformation quantization. In fact the Wigner-Weyl-Moyal 
formulation gives a complete coverage as follows. It is based on the Wigner function (WF) 
which is a quasi probabiity distribution in phase space defined via 

(4.1) f(x, P ) = ^- [ dyr((x-^y]e- i y^(x+^y 



2vr J \ 2" ) V 2 

Here one has (CI) J dpdxf(x,p) = 1 and in the classical limit h — > / reduces to the 
probability density in coordinate space x (usually highly localized) multiplied by delta 
functions in momentum; thus the classical limit is "spiked" . The WF is manifestly real and 
constrained by the Schwartz inequality to be bounded with — {2/K) < f(x,p) < (2/K) (the 
bound disappearing in the spikey classical limit). Projection in x or p leads to marginal 
probability densities, namely, a spacelike shadow f dpf(x,pp(x) or a momentum space 
shadow J dxf(x,p) = cr(p). The WF can and most frequently does become negative in 
some regions of phase space. Nevertheless WF is a distribution function providing the 
integration measure in phase space which yields expectation values from phase space c- 
number functions. Such functions can be classical functions but in general are associated to 
suitably ordered operators via Weyl's correspondence rule. Given such an operator ordered 
via 



ir(p— p)+ia(i~ x) 



(4.2) <S(y,p) = — ^— J dTdadxdpg(x,p)e 

the corresponding phase space function g(x,p) (classical kernel) is obtained via (C2) p — > p 
and y — > x. The operators expectation value is then a phase space average (C3) < >= 
J dxdp f(x,p)g(x,p). The classical kernel is often the unmodified classical expression such 
as H = (p 2 /2m) + V(x) but it contains h when there are ordering ambiguities (see below). 
This operation corresponds to tracing with a density matrix as indicated below. 

The dynamical evolution is specified by Moyal's equation which extends the Liouville 
theorem of classical mechanics (CM), namely dtf + {/, H} = 0, and is given by (C4) dtf = 
(l/ih)[H ★ / — / * H] where (C5) * ~ exp[(ih/2)( d x d p — d p d x ). The right side of 
(C4) is of course the Moyal bracket and is the essentially unique 1-parameter associative 
deformation of the Poisson bracket (cf. |33| F53 ) . In practice evaluation can be expressed 
through Bopp operators in the form 



( I-'!) f(x,p)-kg(x,p) = f [x + ^-d p ,p- y d x ) g(x,p) 



The equivalent Fourier representation of the star product can be expressed via 

(4.4) f* g = —L- [ dudvdwdz f(x + u,p + v)g(x + w,p + z)e {2i/h)( - uz - vw) 
(Jmy J 
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which exhibits noncommutativity and associativity. There is a complete isomorphism be- 
tween star multiplication and operator multiplication indicated in 

(4.5) 21<B = — ^ [ dTdadxdp{a * &) e ir(t.-p)-H«r(r-x) 

(2viT J 

One sees also from (|4.4j) that (C6) j dpdx f * g = J dpdx fg = J dpdx g* f ■ Note that the 
Moyal equation is necessary but does not suffice to specify the WF for a system (e.g. f(H) 
commutes with H). 

Static or stationary WF's obey more powerful stargenvalue equations (cf. |261 1271 13*2*] ) 

(4.6) H(x,p)* f{x,p) = H Ix + — d p ,p - — d x J f[x,p) = f(x,p)*H(x,p) = Ef(x,p) 

where Hip = Eip and this amounts to a complete characterization of the WF's. Indeed, 
using a simple Hamiltonian p 2 /2m + V(x) (without essential loss of generality) one proves 
(cf. [65 ) that for real f(x,p) the Wigner formula (|4.1() for pure stationary eigenstates is 
equivalent to compliance with the stargenvalue equations ("4.6JI along with / ★ H = Ef. 
Conversely the pair of stargenvalue equations for f(x,p) = J dyexp(—ipy)f(x,y) leads to 
/ = (l/2n)ip* {x — (h/2)y)ip(x + (H/2)y). There are also a number of special properties for 
pure state f. Thus from (C7) f*H*g = Eff-kg = E g f*g so for E g ^ Ef one has (C8) f*g = 
0. Moreover f or / = g one gets then (C9) f*H*f = E f f*f = H*f*f so (CIO) /★/ oc / 
and in fact (Cll) f a *fb = (^/fydabfa- Here the normalization is important since it prevents 
superposition (which is handled differently as in the density matrix formulation). Note also 
by virtue of (C6) for different stargenfunctions one has (C12) J dpdx fg = so one must 
go negative to offset positive overlap (a virtue of negativity). Further note that integrating 
()4.6|) yields the expectation of the energy (C13) J H(x,p)f(x,p)dxdp = E J f dxdp = E 
and from (Cll) we get (C14) J f 2 dxdp = l/h. 

Next note (C15) < g* * g >> which leads to the uncertainty principle. Indeed 

(4.7) / dpdx(g**g) =H [ dxdp(g**g)(f*f) = h f dxdp(f*g*)*(g*f) =H f dxdp\g*f\ 2 



To produce Heisenberg's uncertainty principle one chooses (C16) g = a + bx + cp for arbi- 
trary constants a,b,c € C. The resulting positive semi- definite form is then (C17) a* a + 
b*b < x*x > +c*c < p*p > +(a*b + b*a) < x > +(a*c + c*a) <p> +c*b < p-kx > +b*c < 
x-kp >> 0. The eigenvalues of the corresponding matrix are then non-negative and so must 
be the determinent. Some calculation then leads to (C18) AxAp > h/2 (cf. [HE] and see 
I23I for details in calculation). 



Going now to time evolution the Moyal equation (C4) is formally solved by virtue of 
associative combinatoric operations completely analogous to Hilbert space quantum me- 
chanics (QM) through definition of a star-unitary evolution operator (star exponential) in 
the form 

(4.8) U.(x,p, t) = el tH/h = 1 + (it/h)H(x, p) + &^H * H + * H * H + • • • 
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Given the WF at t = the solution to the Moyal equation is then (C19) f(x,p,t) = 
U~ (x,p, t) -k f(x,p, 0) -k U*(x,p, t). For variables x,p this collapses to classical trajectories 

(An\ ■ x-kH-H-kx . p-kH-H-kp 

(4.9) x = = d p H; p = = -d x H 

in in 

Thus for the harmonic oscillator (C20) x(t) = xCos(t) + pSin(t) and p(t) = pCos(t) — 
xSin(t) so the functional form of the WF is preserved along classical phase space trajectories 
via (C21) f{x,p,t) = f(xCos(t) - pSin(t),pCos(t) + xSin(t),0). 

5. MOMENTUM CALCULI 

We summarize first some features of the basic situation. 

EXAMPLE 5.1. Thus go back to Remark 2.1 and the Hamiltonian theory for dKP in the 
form (|2.9j) . written here as (Dl) P = d0. n and X n = —dp0. n . We know KP is equivalent 
to MdKp (Moyal dKP) and examine here the nature of quantizing (Dl) via Moyal- Wigner- 
Weyl (MWW) formulas. In (Dl) the T variable is r n = T' n = nT n where et n ~ T n and 
we recall here the origin of (Dl) from |161 1171 141j . Thus as in Section 2 one arrives at 
d n S = <B„ and we set nH n = 23 n . Then rescaling T n = T' n = nT n and writing d n ~ d/dr n 
now we have (for P = d x S and Q5 n = Q5 n (P,X)) (D2) d n P = d x Q n = dQ n + d P £l n dP 
(in an obvious notation). Then thinking of P = P(X,t) write 

(5.1) P n = 8 n P + dPX n = 8£ln + d P Q n dP + dPX n 

This is then incorporated into a HJ theory with S as a generating function and equa- 
tions (Dl). Now the quantization of the dynamical system via classical MWW meth- 
ods involves Hamiltonians H(X,P) = £l n (X,P) and equations (using r as the evolution 
time and (X,P) ~ (x,p)) (D3) ihX n = {X,£l n } M and (D4) ihP n = {P,£L n } M where 
{f,g}M = f * g — g * f with ★ as in (C5). Hence as in (|4.9|) we have (D5) X n = dp£l n 
and P n = — dxO-n repeated (i.e. the "motion" is along the "classical" trajectories). One 
can ask now what this means in the standard KP theory. Recall from Remark 3.1 how 
the symbol calculus for PSDO is equivalent to Moyal dKP so in some sense we will have a 
Hamiltonian theory in the symbols of KP operators (based perhaps on some sort of relation 
such as (I3~T21 1. ■ 

Now we look at momentum calculi (cf. also [SI 1291 159*] ). The approach of H is based on Lie 
algebras, Poisson structures, and R-matrices; the demands of the general framework adopted 
however seem to restrict severely the range of applicability (e.g. KP is not included). Thus 
we mainly omit this. Another, more flexible, approach is developed in [2*3 (cf. also |59U62j 'l 
One begins with the standard phase space star product 

(5.2) A(x,p) * B(x,p) = e^ d P- d ^A(x,p)B(x,p)\^ =(Xtp) 

with conventional Moyal bracket (D6) {A(x, p), B(x,p)} K = (1/2k)(A*B — B*A). As usual 
one has lim K ^o{A, B} K = {A, B} (Poisson bracket). The star product gives the momentum 
an operator character via (note h/2 ~ k when comparing notations) 

(5.3) p n * p m = p m+n ; p n * f(x) = ]T Q (2 K )"7 (™) * p n ~ m ; 
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n \ n(n — 1) ■ ■ ■ (n — m + 1) f 



1 



jnj ml \0, 

Up to normalization these are precisely the relations satisfied by the derivative operator. 
Let us check here some calculations based on (j.3.6j) where 



_ _ _ _ oo „ s 

(5.4) / *g = fe <d,d p -a p a x ) g = J- « ^(-ly ^^/)^^) 

o s - 

(5.5) p m * <? = E ^pP m ^9 = E ( T J K ^ m "^ (s) ; 

/*f = Ei(- 1 )W = f;( ^ )(-k)V*-Y w ; 

s " ^ ' 

m— s 


Look at some low order terms 1*/ = /, p* f = pf — nf, f *p = pf — nf, and / * 1 = /. 
Then p* f = pf = nf = f * p + 2nf = / * p + (2«)(/' * 1) as in For p 2 we 



have p * f = p > f + i i J K Pf' + K /" witn f *P = P f ~ [ \ J + K f" ■ Hence 
p 2 */ = /V+2« ( J ) = /V+2k ( J ) [/'*?+«/"] = /*P 2 +(2«) ( J ) (/'*?) + 

(2k) 2 /"* 1 = S^o f m ) ( 2K ) m f {m) * P n ~ m ■ Thus seems reasonable and to compare 

with derivation operators note also (k<9/) 2 / = K 2 (f"g + 2/'</ + /</') and (f'K,d)g = fug' 
so 4(k<9) 2 / ~ /(2k<9) 2 + 4k/'(2«9) + {2k) 2 f" . Thus one has 

(5.6) (2«d) n / = E ( £ ) (2K) fc / (fc) (2Kd) n - fe 

o ^ ' 

upon extrapolation. 

Now one defines two classes of Lax operators on the phase space via 

(5.7) L n =p n + ui{x) * p n - x + u 2 {x) * p n ~ 2 + --- + u n (x); 

A n =p n + ui{x) * p n ~ l + ■■■ + u n (x) + u n+1 {x) * p' 1 + ■ ■ ■ 

Thus one has replaced the space of pseudo differential operators by that of polynomials in 
momentum which inherits an operator structure through the star product and defines an 
algebra. This will be called the momentum algebra M n and one notices that this is different 
from the concept of pseudodifferential operators (PSDO) with the coefficients taken from 
the Moyal algebra of [57]. All of the properties of PSDO carry through with suitable 
redefinitions. In particular thinking of the residue as the coefficient of the p _1 term one 
gets (D7) Res{A, B} K = (d x C) exhibiting the residue as a total derivative. Consequently 
one can define (D8) Tr(A) = J dx Res(A) which is unique (with the usual assumptions 
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of asymptotic decrease) and satisfies cyclicity. For a general Lax operator A n one checks 
immediately that 

(5.8) V^M}; W> 

defines a consistent Lax equation provided m = 0, 1, 2 and the projectors are defined with 
respect to the star product (the differencein ordering here can be adjusted via tk — > — tj~ if 
desired). Note here A fc//n = A 1//n * • • • * A 1 /™ with k factors with the n th root determined 
recursively. The projection with m = is denoted by ( )+ and will be referred to as 
the standard Moyal-Lax representation (the others are called nonstandard and are not 
considered here). Note that (D9) lim K ^o(A* A') + = (AA')+ where the factors on the right 
are phase space functions (not operators). Thus one can go the Lax representation of the 
dispersionless limit in a natural manner (cf. |16| I17[ I59|). In fact in this limit one gets 
(D10) dtA n = {A n ,(An //n )+ with the standard Poisson bracket. Further one determines 

conserved charges via (Dll) H k = Trk k J n (k^£n), proves that different flows commute, 
and defines Hamiltonian structures in a straightforward manner. For illustration consider 

EXAMPLE 5.2. In the KdV hierarchy is developed via the Lax operator L = p 2 +u{x) 
where (D12) (L 3 / 2 )+ = p 3 + (3/2)it * p + (3k/2)u^ (uW = d x u - we have changed the 
coefficient of to agree with calculations below in (D25)). One gets then 

(5.9) d t L = |l, (i 3/2 ) + | => d t u = - (ku 3 + j 
The first few conserved quantities are (not checked) 

(5.10) Fi = Tr(LV 2 ) = J dx(u/2); H 2 = Tr(L 3 / 2 ) = J dx(u 2 /4); 

H 3 = Tr(L 5/2 ) = J dx(inu (2) u + n 3 ) 

The commutativity of flows follows directly from the Moyal-Lax representation. ■ 

EXAMPLE 5.3. The conventional Lax equation in standard representation (D13) dt k = 
[(L fc / n ) + , L] (L a PSDO) resembles a Hamiltonian equation with (L fc / n ) + as a Hamiltonian. 
However such a relation cannot be further developed in the language of PSDO. In contrast 
consider the Moyal-Lax representation with an arbitrary flow in the KdV hierarchy described 
by (D14) dtL = {L, (L^ 2n+1 ^ 2 ) + } K (L as in (|5.4[)1. Then consider an action of the form 
(D15) S = J dt{p*x — {L^ 2n+l ^/ 2 ) + ). It is important to remember here that L = L(p, x) but 
does not depend on time explicitly. Thus one can think of (L^ 2n+1 ^ 2 ) + as the Hamiltonian 
on the phase space. That this is valid follows from the Euler-Lagrange equations ((2n + 
l)/2 = a - cf. also Example 4.1) 

(5.11) x = = {x, (L a ) + } K ; p = "^f^ = {P, (L a ) + } K 

Further since L is a function on the phase space one has (D16) dtL = {L, (L^ 2n+1 ^ 2 ) + } K 
so that the Moyal-Lax equation is indeed a Hamiltonian equation with L" playing the role 
of Hamiltonian. This procedure also goes through for the nonstandard representations. ■ 
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The Moyal-Lax representation has now an advantage in that one can go to the disper- 
sionless limit of an integrable system by simply taking the limit k — > and, while the Lax 
representations for various dispersionless integrable models are known (cf. |111I161IT7| ). the 
determination of Hamiltonian structures (at least the second) from such a Lax representa- 
tion has often been open. The Moyal-Lax representation provides a solution to this problem 
in a natural way (cf. [I'll I29| for details). 

REMARK 5.1. An approach similar to (2H1 was developed in [HH] (cf. also |6*2"] ) 
and we sketch a few points here. One considers an algebra of Laurent series of the form 
A = {A; A = X^-oo a iP 1 } w hh coefficients depending on t\ = x, t2, ■ ■ ■ ■ A can be decom- 
posed as A = A>fc A<fc for k = 0, 1, • • ■ where e.g. A>^ = {A = Yli>k a iP l }'i the notation 
A + ~ A>o is used as before. Evidently A is an associative but noncommutative algebra un- 
der the Moyal star product and one defines Res(A) = a_i with trace as Tr(A) = J Res(A). 
There results 

(5.12) / Res(A * B) = j £ J —^^— { a ihj )^ = £ J aA+1 

One notes that this is the same as in the dispersionless limit k — > and because of this the 
Hamiltonian formulation for the Moyal KdV (for example) becomes possible. Indeed using 
(t5~T2l one shows that Tr{A, B} K = and Tr(A*{B, C} K ) = Tr({A, B} K *C). To see this use 
(|5.9jl to replace the Moyal star product within the trace by the ordinary multiplication. Now 
given a functional F(A) = f f(a) one defines a gradient via (D17) cIaF = ^(<J / / 5di)p~ 

where the variational derivative is defined via (D 18) 5f/5ak = iy(d l -(df/da^) where 

af = (d i ■ ofc) and d ~ d/dx. Note d- f = f = df/dx and df = fd + f. The Moyal KdV 
hierarchy is defined by the Lax equations 

(5.13) ^- = {(L l ^) k + , L} K = {L, } k 

where {L l l n ) k + = L l l n * ■■■* L l / n (k times) and here L = p n + Eo _1 ^ with Ll,n = 
P + ^2 a iP~ l is the n th root (D19) L = L l / n * ■ ■ ■ * L x l n (n times). By definition of the 
Moyal bracket the highest order in p on the right side of H5.1UJI is n — 2 so u n _i is trivial 
in evolution equations and can be dropped in the Lax formulation (this changes in the 
Hamiltonian formulation however). The simplest example is written out in [^2], namely 
n=2 with L = p 2 + u and L 1 / 2 = p + J2T a iP~* and ( D2 0) ai = tt/2; a 3 = -u 2 /8; a 5 = 
u 3 /16 + (k 2 /8)(u 2 — 2uu x ); ■ ■ ■ and the first few Lax flows are (D21) Ut 1 = u x ; Ut 3 = 
(3/2)uu x + k 2 u xxx ; ■ ■ ■ . This set of equations forms a Moyal KdV hierarchy which can also 
be obtained from reduction of the Moyal KP hierarchy (cf. |11[ I12j or noncommutative zero 
curvature equations (cf. 42 J. When k = all higher order derivative terms disappear and 
the Moyal KdV reduces to the dispersionless KdV hierarchy (cf. JJ)- ^ n this sense the 
Moyal parameter k characterizes the dispersion effect. On the other hand when k = 1/2 
the Moyal KdV hierarchy is the ordinary KdV hierarchy and this is due to an isomorphism 
of Moyal KP to ordinary KP at k = 1/2 (cf. |11 [ 1361 155] and Section 3). ■ 

We concentrate now on Examples 4.1 and 4.3 and for simplicity begin with KdV. The 
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theme is somehow to exhibit the quantum mechanical nature of KdV via the dynamical 
system 1)5.1 Jl or 1)5.11)1 . The momentum algebra as in this section which leads to 1)5.11)1 and 
(D16) is very attractive and we spell out more details. We go first to [55] and extract 
some specific formulas to enable the study of (D22) dkL = {L^ 2n+1 ^ 2 , L} K (tk —>■ —tk in 
(D16)). We will not emphasize the bihamiltonian structure. Then for L = p 2 + u one has 
(D23) L 1 / 2 = p + J2T a iV~ { with ( cf - ( D12 )) 



(5.14) L 3/2 = L-k L 1/2 = (p 2 + u) * (p + <HV l ) = p 3 + p 2 * a^ 1 + p 2 * a 2 p~ 2 + ■■■ + 

l 

+u -kp + u-k a\p~ l + u-k a 2 p~ 2 + ■ • • 
From we have a 2 k = and (cf. (D20)) 

(5.15) ai = -u; a 3 = ~^u 2 ; a 5 = -^u 3 + ^k 2 (uI - 2uu xx ); 

5 5 1 

«7 = ~Y28 n4 ~ i6 k2 ^ MU ^ ~ 2u 2 v>xx) ~ g« 2 ("L ~ 2u x u xxx + 2nn (4) ); • • • 



(5.16) u tl = u x ; u t3 = ^uu x + k 2 u xxx ; u t5 = —u 2 u x + -k 2 (uu xxx + 2u x u xx ) + k 4 u (5) - 



The first few Lax flows are given by 

^uu x + k 2 u xxx ; u t5 = ^u 2 u x + ^/ 

The comparison here is with (D12) (as changed from 29 ). Thus the only plus terms 
involve the composition (since a 2 = 0) 

(5.17) p 2 -k a\p~ l = a\p + 2^d x ai + K 2 p~ 1 d x a\ 

so (D24) L 3 / 2 = p 3 + aip + 2nd x a\ + u*p = p 3 + (l/2)up + ku x + u-k p. We know that 

u-k p = pu — ku x from calculations after 1)5.5)1 so (D25) Z/^ 2 = p 3 + (1/2) [ti *p + ku x ] + 
Kti^: + n -kp = p 3 + (3/2)u *p + (3/2)ku x . 

REMARK 5.2. We see via Examples 4.1 and 4.3 that KP or KdV for example have a 
distinctly quantum mechanical flavor via the phase space dynamics (|5. 11|> (note again that 
it doesn't matter whether we take {L°L , L} or {L, L"} to express the time variation since it 
is simply a matter of time reversal) . Thinking of KdV for simplicity, as a result one say that 

all du/dt 2n +i aarise from QM equations (D26) d 2n+ \ = {L, L^ 2n+1 ^ 2 }. Further, given the 
dependence of the ai in (D20) or 1)5. 15|) on u one has a direct computation for da 2 i + \ / dt 2n+ \ 
(this also results from (D26) upon writing out terms). The QM nature of the flows arises 
explicitly via the presence of k in the equations and this reduces to "classical" behavior 
when k — > (where classical here is unrealistic since it ignores dispersion however). Thus 
dispersion for example appears as a quantum phenomenon; QM represents a smoothing or 
calming factor which eliminates some wave caustics, breaking, etc. Since one knows also 
that q ~ exp(n) is characteristic of the emergence of q in various q-versions of quantum 
phenomena (cf . ^2] ) we have here another physical motivation for q which makes the study 
of q-QM seem more meaningful (one can ask here also about whether a k quantization as in 
Moyal is just a first order version of a q = exp(n) discrete quantum theory). On the other 
hand water waves for example are eminently classical macro-phenomena so one can ask why 
they should obey QM rules of behavior. As an aside we mention, following |29| . that upon 



14 



ROBERT CARROLL UNIVERSITY OF ILLINOIS, URBANA, IL 61801 



constructing the bihamiltonian theory associated with (D26) one finds that k is directly 
related to conformal field theory (CFT) and to the central charge of the second Hamiltonian 
structure of KdV (which in turn is related to a standard Virasoro algebra - not q-Virasoro 
- cf. HUHUE1). If one thinks of KdV (or KP) then as a QM extension of dKdV (or dKP) 
then some real world macro-phenomena are essentially quantum mechanical in nature. This 
is perhaps a viewpoint to be further explored both philosophically and technically. Note 
that (elementary) particles have only one QM time evolution to drive them but in some 
sense "fluids" seem to have many, whose nature could be further examined with the study 
of higher order evolutions in KdV for example. ■ 

EXAMPLE 5.4. Following |291I59| we compute the first terms of the "quantum" equation 

d t L = {L 3 / 2 ,L} K with L 3 ' 2 = p 3 + {3/2)ukp + {3/2)ku' (cf. (D12) and (D25)). Thus first 

we get (since u* v! = uv!) 

3 3 3 3 

(5.18) T = (p 3 + -u-kp + 2 K-U ') * (P 2 + u ) - (P 2 + u)-k{p 3 + -u*p + 2 Kn ') = 

^ ^ S ^ *^ / 2 r ^ 2 / \ *^2 / S i 

= P *U + — U-kp + —U*p-kU + —KU -kp — [—p * (U * p) + — p kr KU + U -k p + —U * U -k p\ 

Note now from (|5.5j) that (D27) p 3 *u = k 3 u'" + 3n 2 pu" + 3np 2 u' +p 3 u while (D28) u-kp 3 = 
—K 3 u"'-\-3K 2 p—3Kp 2 u'+p 3 . Hence (D29) p 3 -ku—u-kp 3 = 2k 3 u"' +6np 2 u' . Next consider H = 
(3/2)(u*p 3 — p 2 -ku-kp) with u-kp = p*u — 2ku' so (D30) S = (3/2)(u-kp 3 — p 3 *u+2Kp 2 -ku f ) = 
-(3/2) (2k 3 -u'" + 6kp 2 u') + (3/2)2np 2 *u' = -9Kp 2 u' + 3n(p 2 u' + 2Kpu") = -%np 2 v! + 6npu" . 
Finally (D31) (3 /2)(u* p-k u — u* u* p) = (3/2) (if ★ (u-kp + 2kv!) — u-ku-kp) = 3niu-ku and 
(D32) (3/2)k(u'*p 2 -p 2 *u') = -6npu". Combining then D27) - (D32) we get (D33) T = 
2k 3 u"' + 3nu-ku' so the Lax equation then decrees that (D34) Ut = (1/2k)(2k 2 u'" + 3kuu') = 
{3/2)uu ! + kv!" which is the standard KdV form with the dispersion coming from k. ■ 

6. REMARKS ON Q-THEORIES 

QKP (and qKdV) can be developed in a hierarchy form following [U IH EH ESI 133 EH 
I4fl[ I60| I61j . This is pursued in connection with Hirota formulas at some length in |19j and 
we remark here only that the resulting qKP or qKdV equations are very complicated due 
to formulas of the form 

(6 .i) „ = n = (i - d) (wm-QiM _ (Oj^V + m d ( b jh) + 



instead of the classical u = d 2 log(r). Similarly for qKdV one has a difficult formula 
(A31) u = d q dilogr(x,t)DT(x,t). 

REMARK 6.1. In this direction if one actually writes out a qKdV equation for example 
from the hierarchy picture it will have the form (cf. |19j ) 

(6.2) d t u = (d 3 u) + w 2 (d 2 u) + wi{d q u) - [(d 2 w ) + ui(d q w )] 
where (setting (El) u\ = (q — l)xu = (1 + D)sq and u = s\ + Ds\ + Sq + d q so) 

(6.3) w 2 = D 2 s + u\ = D 2 s + Ds + s ; 

w 1 = {q+ l)(Dd q s ) + D 2 Sl + [(Ds ) + s ](Ds ) + u; 
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w = dgS + (q + T)(Dd q si) + uid q s + ui(Dsi) + us + D 2 s 2 

(cf. |12H2L)llHn| for other forms of q-equations or noncommutative integrable equations). It 
appears therefore that after expressing e.g. so = (1 + -D) _1 (g — l)xu with (E2) (1 + D)^ 1 ~ 
X^o°( — D) n formally the qKdV equation will have an infinite number of terms. A similar 
kind of equation for "qKdV" was derived in |21| (cf. also [21]), by use of a version of q- 
Virasoro, in the form (E3) ut + cd q (D + D~ 1 )~ 1 d q u + d q (uDu) + D~ l ud q D^ 1 u = where 
here d q f = [f(qx) - f{q~ l x)]/{q - q^)x. ■ 

One recalls here that d q f(z) = [f(qz)-f(z)]/(q-l)z with (E4) 8 q (fg) = d q {f)g+ T {f)d q g 
where r(/)(z) = f{qz) (we use r and D interchangably now and note that d q T = qrd q ). 
PSDO are defined via an equation (E5) A(x, d q ) = X^-oo u i( x )dq with d q u = (d q u) + r(u)d q 
and one has (D q ~ d q ) 



(6.4) a- 



Recall also 



k>0 



-fc(fc+l)/2 (r -fc-l ( ^ u)) ^-fc-l. 



E 

fe>0 



n 
k 



{r n - k {d k q u))d n q 



(6.5) 



9-1 



m 
A- 



(m) q (m — l) q ---(m — k + l) q 



(l) q (2) q ---(k\ 



d q + ui(z) + u 2 (z)d q 1 



+ 



There are then q-analogues of the Leibnitz rule etc. and for L q 
uz(z)d~ 2 + • • • and one has q-KP via (8L q /dt m ) = [(L 1 q n ) + , L q ] where u\{z) has a nontrivial 
evolution. In accord with Section 3 we should now represent the ring or algebra 2t g of qPSDO 
symbols via a product as in say (B7) and thence provide expressions for deformation thereof. 
The X and P variables should come from the phase space for dKP. Evidently the qPSDO 
symbols will involve a variation on (|3.21|) and one can utilize techniques of pil3lHi rL2"ll5L)lHU] 
(cf. also |13M19ll2"U] for q-formulas) for calculations. In this direction one finds (heuristically) 

-> PSDO and % q correspond symbolically via d ~ 
^ r factors along with q-subscripts; in particular 
is needed. When commutators are also envisioned 



that the ring or algebra calculi of 21 
D q = d q and suitable insertion of D 



(E6) (d^e 



(i-k + l) q (rO % 



have different degrees, etc. but one 



individual terms may differ because e.g. brackets 
notes that e.g. 

(6.6) [ad, bd] = (ab' — ba')d; [ad q , bd q ] = (ab' q — ba' q )d q + (arb — bra)d q ; 

(ab' — ba')d — > (ad q b — bd q a)d q = [a(b' g + rbd q ) — b{a' q + rad q )]d q 



REMARK 6.2. Given L q = D + q + b (t) + J2T h -t D q i for q KP and L = d + 
^2^° f3-id~ l for KP, with nontrivial evolution of bo one cannot perhaps expect KP and qKP 
to be "isomorphic" but as indicated above they do correspond symbolically so there seems 
to be no reason to regard one as more "quantum" than the other (except perhaps the fact 
that the q-theory represents a discretized version of the other and such discretization may 
be an essential quantum signature - not yet authenticated in conventional QM such as the 
Moyal quantization). There is also another way to view qKP due to E. Frenkel (cf. |35j ) 
via (E7) qKP : L q = D + a (t) + a_i(t)Z) -< where Df(t) = f(qt) (recall here that 
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{q — l)xD q f = (D — l)f). It is shown in e.g. [21E1 (cf. also |12|) that there is an isomorphism 
mapping qKP or qKP into the discrete KP hierarchy which exhibits their equivalence via 

" k + i ~ 

( 6 - 7 ) a ^ = E ( _ v( *i) a i)k b >°+i(y) 

0<k<n-i y yVH ' H ' 

We refer to [21 [H] for details on discrete KP which in fact is equivalent to the 1-Toda lattice 
(cf. [21 121)- The correspondence to qKP can be seen best through the notation of [SS] where 
one works with Toda Lax operators L = A + ^u n+ i(e,f,[,s)A" n with A ~ exp{ed s ) 
(cf. also We note also connections of discrete KP to algebraic equations over finite 

fields in [3] and an interesting paper [3J on formulating quantum mechanics with difference 
operators. ■ 
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